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Hybrid quantum systems usually consist of two or more subsystems, which may take the advantages of the 
different systems. Recently, the hybrid system consisting of circuit electromechanical subsystems have 
attracted great attention due to its advanced fabrication and scalable integrated photonic circuit techniques. 
Here, we propose a scheme for high fidelity quantum state transfer between a superconducting qubit and a 
nitrogen-vacancy center in diamond, which are coupled to a superconducting transmission-line resonator 
with coupling strength gi and a nanomechanical resonator with couphng strength g2, respectively. 
Meanwhile, the two resonators are parametrically coupled with coupling strength /. The system dynamics, 
including the decoherence effects, is numerical investigated. It is found that both the small (/<S; {gi ,g2 } ) and 
large ij^lgi^gi}) coupling regimes of this hybrid system can not support high fidelity quantum state 
transfer before significant technique advances. However, in the intermediate coupling regime (/~ gi ~ ga); 
in contrast to a conventional wisdom, high fidelity quantum information transfer can be implemented, 
providing a promising route towards high fidelity quantum state transfer in similar coupled resonators 
systems. 

A rapid processor and reliable memory are indispensable components for building quantum computers. 
Recent progress in mesoscopic objects, e.g., superconducting qubits, shows that this kind of artificial 
systems can couple strongly to electromagnetic field and enable fast quantum logic gates'. However, they 
possess relatively short coherence time. On the other hand, microscopic systems, e.g., nitrogen-vacancy (NV) 
centers in diamond", naturally have rather long coherence time. Therefore, it is attractive to combine the two types 
of systems by employing their hybrid architecture''. In order to take the advantages of both, quantum state transfer 
(QST) between the two sides is crucial. To this end, some strategies have been proposed with hybrid systems* '", 
where a superconducting qubit acts as a processing unit and the other ones serve as memory units. 

There are generally two kinds of hybridity: coupling the two elements via a quantum bus' or directly*. First, a 
superconducting qubit can be coherently coupled to an ensemble of NV centers in a diamond via a super- 
conducting transmission-line resonator, which acts as a quantum bus'"'. However, the NV center ensemble 
usually has relatively much shorter coherence time comparing with single NV center and the coupling of single 
NV center to the transmission-line resonator is too weak to be usable, usually under 0.1 kHz'. Meanwhile, a 
superconducting flux qubit can be directly coupled to an NV center with a coupling strength on the order of 
10 kHZ" '". But, such schemes are difficult to extend to couple long distance qubits as the superconducting qubit 
can not be too large. Therefore, to be potentially extended to quantum networks", many QST schemes have been 
investigated in systems consisting of two coupled resonators'^"^'. Typically, an NV center and a superconducting 
qubit can be coupled to a nanomechanical resonator^"'"^'* and a transmission-line resonator', respectively. On the 
other hand, the two different resonators can be parametrically coupled to form an electromechanical architec- 
ture^^. With rapid experimental progress^^"^", this typical system can have many potential applications for hybrid 
quantum information processing" '" and entanglement^'^ 

In this paper, we propose a scheme to realize high fidelity QST between a superconducting qubit and an NV 
center, which provides the possibility of storing a state of a superconducting qubit into the spin-based quantum 
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memory via the electromechanical system. Previous works in similar 
coupled resonator systems, such as Refs. 17-21, were mainly focusing 
on the large inter-resonator coupling regime (we wUl specify differ- 
ent regimes later), where analytical result can be obtained^'. 
However, after taking the decoherence effects into consideration, 
we investigate the QST dynamics numerically and find that both 
the small and large coupling regimes of this system can not support 
high fidelity QST with current technology. The main result of present 
work is the comprehensive study of QST in the intermediate coup- 
ling regime, which is lacking in previous studies'^"^'. Qur motivation 
also relies on the fact that the intermediate coupling regime can be 
realized more easily as it requires smaller inter- resonator coupling 
comparing with large coupling regime, and thus our result is useful 
and can be tested in a near future experiment. Interestingly, the main 
finding here is that fast and high fidelity QST is possible in the 
intermediate coupling regime, even better than that of the large 
coupling regime. The physics behind this surprising result is that 
the decoherene has less effect in the intermediate coupling regime 
since the time of QST is shorter than that of both the small and large 
coupling regimes. Therefore, our result provides a promising route 
towards high fidelity QST in similar systems with coupled 
resonators. 

Results 

The hybrid system. The hybrid system consists of a superconducting 
qubit (NV center) coupled to a transmission-line (nanomechanical) 
resonator and the two resonators are parametrically coupled. Our 
final goal is to transfer the quantum state of a superconducting qubit, 
the processor, to the NV center, the quantum memory, via the two 
coupled resonators. Since the QST process interchanges the informa- 
tion between the two subsystems, the main results we obtained in the 
following should also be valid for the reverse transfer, i.e., transfer 
the state of the NV center to the superconducting qubit. For 
convenience, hereafter, the subscript 1(2) stands for the supercon- 
ducting qubit (NV center) and a{b) for the transmission-line 
(nanomechanical) resonator, respectively. A superconducting qubit 
can be treated as a two-level system and the transmission-line 
resonator can be modeled as a single Bosonic mode^^"": = 
w^a^a, where we assume H = 1 hereafter, to^ and a (a^) are the 
mode frequency and the annihilation (creation) operator of the 
resonator, respectively. The superconducting qubit can be 
capacitively coupled to the transmission-line resonator and the 
interaction HamUtonian can be written as^'^ '^ 



(1) 



where c* = (cf ± ia'')/! are the raising/lowering operators of the 
superconducting qubit and the coupling strength can be tuned by 
an external magnetic flux of the qubit loop^^'"". 

Meanwhile, the ground state of a negatively charged NV center is 
spin triplet = 0 and = ±1, and there is a zero-field splitting 
between them due to spin-spin interaction. The two states | ± 1) are 
coupled to the |0) state by different polarized driven fields with the 
same detuning A > 0 and strength Q'"'. In a frame rotating with the 
driven frequencies, the Hamiltonian of the NV center can be written 



Hnvc = 



A|i)(/| + -fi(|0)(i| + |/)(0|) 



which couples |0) to the bright state \b) = + \l)) ^ ^2, while 

leaves the dark state = — 1) — |l)) ^ \fl uncoupled. The eigen- 
states of the NV center are |m) = cos i9jO) — sin i9|fc), and je) = 
cos i9|fo) + sin i9|0) with tan(2,9) = \fl£L / A and the corresponding 



eigenvalues are labeled as coj^^.e- The Hamiltonian of a nanomecha- 
nical resonator can be written as Hh = tO[,b^b with co^ and b (fo"'") being 
the mode frequency and the annihilation (creation) operator of the 
nanomechanical resonator, respectively. A magnetic tip is attached 
to the nanomechanical resonator at a distance h ~ 25 nm above a 
NV center, and thus creates a strong coupling between the two ele- 
ments. The mechanical resonator can be a doubly clamped beam or a 
cantilever, and the magnetic tip can be attached at the position of an 
antinode and the undamped end, respectively. The superconducting 
qubit is far away from the tip (comparing with h), and thus feels 
negligible magnetic field. When aj(, = oj^^^ = Wd ~ C0m> the inter- 
action between the NV center and the nanomechanical resonator 
reads^^ 



Hb,2=g2(bal + b\, ], 



(2) 



where 02 =\m){d\, (^2~{'^2)'' ^^'^ ^2("CW(,) is the coupling 
strength. 

Furthermore, we consider a capacitive coupling between the 
microwave transmission-line and the nanomechanical resonator^^'^". 
The mechanical resonator is driven by an ac voltage, which matches 
the energy difference between the electrical and mechanical systems, 
and thus couples the two. Then, an effective linear coupling between 
the two resonators can be induced, and the effective electromech- 
anical interaction Hamiltonian can be written as 



(3) 



where the coupling strength / determines the speed of the energy 
transfer between the two resonators. The nanomechanical resonator 
here can be a doubly clamped beam^** or a metallic membrane^\ 
However, the direct capacitive coupling strength of the microwave 
resonator and the mechanical beam is usually very small in typical 
experiment^**. Fortunately, one can replace the coupling capacitor by 
a superconducting qubit, which can greatly enhance the coupling 
strength'*". Alternatively, one can use a metallic membrane to serve 
as one of the two parallel metal plates of the coupling capacitor*'*', 
where one can readily obtain / ~ g2. In this case, the coupling 
between the membrane and the NV center can be induced by strain'*", 
where strong coupling can be achieved. Therefore, the total 
Hamiltonian for the hybrid system is given as 



(4) 



The competition of different couplings of this Hamiltonian will yield 
three coupling regimes. As /determines the time scale of the system- 
atic dynamics, we will call /> {^1 ,^2 } . {gi ,^2 } , and 7 ~ ^1 ~ ^2 as 
large, small, and intermediate (inter-resonator) coupling regime, 
respectively. In addition, to get coherent dynamics, the coupled sys- 
tems described in Hamiltonian in Eqs. (1) and (2) need also to be 
operated in the strong coupling regime. 

The quantum state transfer. We now study the QST from a 
superconducting qubit to a NV center. We consider the zero- and 
one-excitation subspaces, i.e., spanned by the basis of {|0m)i2|10)ai„ 
|0m)i2|0l)ai„ I I»i)i2i00)a(„ |0d)i2|00)ai,}. The basis vectors and their 
corresponding populations are labeled as \(^„) and P„(n = 1,2, 3, 4), 
respectively. For an initial state ji//), = |i/'i)|m)2|00)„f, with ji/'j) = 
cos 0\ 1)1 -I- sin 0|0)i, the goal of the QST is to obtain a final state of 
|i/')/= |0)i(cosO|d)2 + sin^)|m)2)|00)„i,- Note that |0)i|m)2|00)„i, will 
not evolve under the act of Hamiltonian (4) and the QST completes 
when the excitation initially in the superconducting qubit (P3 =1) 
transfers to the NV center (P4). We estimate the performance of our 
scheme by the conditional fidelity defined by F = {i/'i|pa|i/'i), with Pa 
being the reduced density matrix of the NV center from the final 
state, under decoherence rates c using the Lindblad master equation, 
see the method section for details. 
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Figure 1 | Fidelity as a function of the reduced time of tl (Itt) in the small 
coupling regime with / = 0. 1 for different ^, other parameters are ff = nli, 
C= 0.001, and gi = g2=l. 

The small coupling regime. In this regime, the dynamics can be 
understood as two weak coupled subsystems with each having a 
Jaynes-Cummings type interaction. The dynamics of the system is 
characterized by a fast oscillation behavior with frequency g=gi =g2 
accompanied by a slow oscillation envelope with frequency J/2. 
Choosing 6 = n/4, the fidelities are plotted in Fig. 1 with different 
^. This transfer is completed at instants T„ = {2k + l)n/J with k 
being a nonnegative integer. However, the maximum fidelity with ^ 
= 0.03 is less than 70%. The relatively large infidelity is mainly due to 
the decoherence effect since the transfer time is very long because of 
the very small /. To obtain a fidelity about 95%, one needs to reduce 
the decay rates to C/10, which is very challenging. Moreover, it is 
worth to note that if the symmetric coupling situation gi = g2 is not 
satisfied, the deviation brings detrimental influence to the QST: the 
population maximum wiU decrease rapidly with the increase of the 
deviation denoted by \g2 — gi\/g2- Furthermore, the fidelity is too 
small, which means it is impossible for QST with asymmetric 
coupling in this regime. 

The large coupling regime. In this regime, the system can be 
considered as two subsystems connected by large photon hopping. 
The relatively large coupling between the two resonators services as a 
fast channel to transfer the energy residing in the two qubits. The 
system dynamics can be characterized by two distinct oscillating 
frequencies, fast oscillation with frequency / and slow oscillation 
with frequency^' g^/J. QST can be completed at instants = (fc + 
l/2)nj/g^ with k being a nonnegative integer. The performance of this 
regime is better than that of the small coupling regime, as shown in 
Fig. 2(a). The maximum fidelity of the transfer process is about 75% 
for present experimental parameter of c = 0.03. To obtain a fidelity 
above 90% and 95%, one needs to reduce the decay of the system to 



f/3 and c/10, respectively. In this regime, as shown in Fig. 2(b), 
asymmetric coupling strength will introduce little influence, and 
thus it is more robust against coupling strength deviation than that 
of the small coupling regime. With the increase (decrease) of from 
the symmetric value of 1, the transfer time will be shortened 
(prolonged) and the fidelity of the QST will be gradually increased 
(decreased). 

The intermediate coupling regime. In the following, we will focus 
on the QST performance in the intermediate coupling regime, which 
is the best choice for QST purpose. First, the QST in the small 
coupling regime needs very long time (r„,~ 1//^1/^), and thus 
decoherence will introduce huge errors. Second, the interaction 
between the resonators is usually weak, or at least is the same 
order as the qubit-resonator interaction g. Though we may tune g 
to match the condition of /^>^, i.e., the large coupling regime, the 
expense is that the effective coupling will be decreased and the time 
to complete the QST will also be relatively long: ~ 1 / /eff = 
J/g^^l/g- This means that the time needed for QST in the two 
regimes are much longer than their respective time scale of 1/^, 
and thus decoherence wiU cause intolerable influence in both 
regimes. However, for the intermediate coupling regime / ~ ^, as 
the time needed for QST wiU be T, ~ 1// ~ 1/^, we expect higher 
fidelity QST than that in the two other regimes. 

For verification purpose, we plot the fidelity of the QST process in 
Fig. 3(a) with ^2 =gi = '^>^ = 0.03, and 0.01 < / < 3. It is clear that 
the fidelity presents regular distribution in the whole region and each 
red spot corresponds to a maximum. When / increases gradually 
from 0.01 to 1.2, the red spot moves closer and closer to the vertical 
axis. When /~ 1.2, the fidelity reaches a maximum value with the 
shortest QST time. If one continues to increase / towards / > 1.2, the 
red spot moves away from the vertical axis. It means that the time 
needed for QST will increase gradually and at the same time the 
fidelity wiU decease gradually. Comparing with that in the Fig. 1 
and Fig. 2, the QST here completes in a much shorter time (the 
maximum value of the fidelity here move towards the vertical axis), 
which leads to the increase of QST efficiency and fidelity. Therefore, 
our numerical results show that the QST performance in the inter- 
mediate regime can be greatly enhanced comparing to that in the 
other two regimes. In order to have an insight into the detail popu- 
lation changes, we focus on the red spot closest to the vertical axis, the 
populations for different states are plotted in Fig. 3(b). The amplitude 
of P4 reaches the largest magnitude, implying the initial state fully 
transferred from the superconducting qubit to the NV center. 
Therefore, J ~ g makes the quantum bus to be a expressway and 
the excitation can pass through faster and with higher fidelity. 

Furthermore, we investigate the influence of asymmetric coupling 
gi ^ ^1 to the QST in this regime. Figure 3(c) shows the time- 




Figure 2 I FideUty as a function of the reduced time of t/(2;r) in the large coupling regime with /= 10, other parameters are 0= jzH and^= 0.001. (a) 

For gi = g2 = 1 with different t. (b) For g2 = 1 and J = 0.03 with different gi. 
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dependent population of the target state P4 with C = 0.03, ^2 = 1 > 7 = 
1.2 and 0.01 < ^1 < 3. One can easily find that the amplitude of P4 
enhances and rapidly moves to the left when g2 increases gradually 
from 0.01. This improvement does not means the success of faithful 
QST, due to the additional influence by vertical adjacent peaks. 
However, within the small region of ^1 = 1.1 ± 0.1, the peaks are 
relatively far apart from each other in vertical direction, which makes 
the higher fidelity QST to be realizable. We have confirmed that one 
can perfectly realize the QST with ^1 = ^2 = 1 ^nd / ~ 1 . 1 6, where the 
QST needs less time (which is about 1/10 and 1/5 of the time in small 
and large coupling regimes, respectively), and the QST has the high- 
est fidelity of 96%. We further investigate the 0 dependence of the 
maximum fidelity, as shown in fig. 3(d), which indicates that the 
larger 6 will lead to larger fidelity. 

In conclusion, we have demonstrated that a high fidelity QST 
between a superconducting qubit and an NV center in diamond in 
a hybrid electromechanical system can be realized in the intermedi- 
ate coupling region. We have numerically simulated the evolution of 
whole system and find that QST can be realized much faster and 
more reliable in the intermediate coupling regime than that in the 
small or large coupling regimes. As the above conclusion is model 
independent, our result may be directly extended to similar systems 
consisting of coupled resonators and thus provides a promising way 
for realizing high fidelity QST. 

Methods 

Modeling of decoherence effects. Inevitably, the QST process will suffer from 
decoherence. We simulate the performance of our scheme under realistic conditions 
by considering the decays of the two resonators {Ka and k;,), the relaxation {y'l and 72) 
and dephasing (r^ and of the two qubits. Under these conditions, the whole 
system can be described by the Lindblad master equation 

;=1.2 

where C{A) — 2ApA^' —A''Ap — pA^'A is the Lindblad operator. We next justify our 
chosen parameters. The superconducting qubit in a transmission-hne resonator can 
be protected by the resonator, for a planar transmon qubit, relaxation and coherence 
times of 44 and 20 /(s are reported for a ID resonator^^, which leads to Fi/ln — 8 kHz 
and yjln = 3.5 kHz. The decay rate of the ID transmission-line resonator is around 
Ka/ln'^S.S kHz {lifetime t„ — 45 j^is^^). The coupling between the superconducting 
qubit and the transmission-line resonator is very strong and can be controlled very 



well, and thus we assume can be tuned to arbitrary value we want, which can be 
achieved by using a flux-biased rf SQUID to couple the superconducting qubit to the 
cavity^^. For a temperature of T — 100 mK and the quality factor Q — 10^, the decay 
rate of the nanomechanical resonator is /vf,/27r^2 kHz^^, andgz/ln::^ 115 kHz^^. As 
the decay rates are in the same order Fi ~ )'i ~ k^, ~ k^, we treat them as a same 
parameter C — i<a- Hereafter, all the parameters will be normalized in unit of g2, and 
thus c — 0.03. The energy relaxation and coherence times of the NV center can be 
very long comparing with others^*, and thus ^2 ^^nd r2 is very small. Therefore, we 
have taken these two parameters as the same in the numerical simulations, i.e., r2 — 
)'2 = t, — 0.001, which corresponds to a coherence time of only 1.4 ms. 
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